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A large (infinitely-dimensional) class of completely inte- 
grable (possibly non-autonomous) spin chains is discovered 
associated to an infinite-dimensional Lie Algebra of infinite 
rank. The complete set of integrals of motion is constructed 
explicitly, as well as their eigenstates and spectra. As an ex- 
ample we outline kicked Ising model: Ising chain periodically 
kicked with transversal magnetic field. 
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In the past three decades intricate algebraic techniques 
(under the names quantum inverse scattering or algebraic 
Bethe Ansatz) have been developed 0] in order to con- 
struct Integrable Quantum Many-body (IQM) dynamical 
systems and the associated complete sets of integrals of 
motion. Integrability of a quantum many-body dynami- 
cal systems is defined in a generalized Liouvilean sense; 
namely by the existence of an infinite set of (indepen- 
dent and local) conservation laws. All the so far dis- 
covered IQM systems are 1-dimensional, typically SU(2) 
spin chains or related systems. Quantum integrability 
is nongeneric but, however, of great importance since it 
has been shown recently that the existence of nontriv- 
ial conservation laws generically leads to ideal transport 
properties with infinite Kubo transport coefficients, and 
deviation from quantum ergodicity in general. 

In this Letter we present a new and elementary ap- 
proach to the construction of IQM 1-dim lattice systems. 
It is based on the particular infinite dimensional dynam- 
ical Lie algebra (DLA) generated and represented by the 
essential dynamical observables (in our case it is gener- 
ated by the Ising Hamiltonian ctJctJ_|_]^ and the inter- 
action with the transversal external field cr|) and for 
which the 'transfer matrix' can be explicitly constructed 
from the commutativity condition. We show that any el- 
ement H of DLA may be considered as a Hamiltonian of 
IQM system and construct an analytic DLA valued func- 
tion r(A) of possibly vectorial spectral parameter A £ 
(for some N > 1), comuting with H, [H, T{X)] = 0. T(A) 
is a formal analog of logarithm of the transfer matrix. The 
integrals of motion (conserved charges and currents) are 
derived as coefficients of Taylor expansion of T(A) around 
A = 0. Therefore we have an infinite-dimensional class 
of IQM Hamiltonian systems. Furthermore, real DLA 
of self-adjoint observables generates infinite dimensional 
unitary dynamical Lie group of even larger class of inte- 
grable quantum many-body propagators of possibly non- 
Hamiltonian (non-autonomous, e.g. periodically kicked) 



IQM systems. As an example we work out kicked 1-dini 
Ising chain periodically kicked with transversal external 
field. Moreover, we explicitly calculate the complete set 
of eigenstates and spectra of the conserved charges (in- 
cluding the Hamiltonian). 

We start with Lie algebra il over an infinite spin chain 
spanned by the spatially homogeneous local observables 



(1) 



where aj , s Cz {I = x,2 — y,3 — z}, are the spin variables 
(Pauli matrices) of spin j obeying the standard commuta- 
tion relations [a^ , al] = 26jk(J^a-l = 2i6jk J2s ^prs<Jj, and 



as the maximal number of digits p of some observable 
(|l|) in the expansion of A in terms of basis (Q). We are 
interested in nontrivial infinite dimensional subalgebras 
of il for which the number of elements with order smaller 
than p grows algebraically (as a function of p) and not 
exponentially ('^ 4^) as for if Q. Indeed we found sub- 
algebra 6, which we call dynamical Lie algebra (DLA) 
(essentially generated by Z[3] and ^[n]), and spanned by 
two infinite sequences of selfadjoint observables C/„ and 
14, 



Vn 



^[1(3- 



41], 



^[1(3—1)2], 
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n > 1, 
n = 0, 
n < -1, 

n > 1, 
n = 0, 
n < -1, 



(2) 



for — 00 < n < oo ((3") indicates digit 3 being repeated n 
times), which satisfy the following commutation relations 



[f^m,f^n] — 2z(V4 — n l^i — m), 
[V^,Vn] =0, 



(3) 



The order of observables C/„ and Vn is |n| -)- 1. The cov- 
ering algebra il is equiped with the Euclidean metric as- 
sociated to the bilinear form (scalar product) 
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iA\B)=l\in —trUA^B), 



L2^ 



(4) 



( tr^ is a trace for a finite system of size L) with respect 
to which (|l|) is an ortho-normal (ON) basis. Further, C/„ 
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and Vn form ON basis of DLA 6 in the same metric. 
Note that (Q) is invariant with respect to adjoint map, 
iadA)B = [A,B], namely ((adAt)B|C) = iB\{adA)C). 

Conservation laws in general autonomous case: Let us 
assume that the Hamiltonian H and the logarithm of the 
transfer matrix T belong to DLA 6. We write 



TV 
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T= ^ {amUm + bmVm)- 



(5) 
(6) 
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where the Hamiltonian H has a finite order M :— 
max{m+,m_}. The commutativity relation [H,T] — 
gives the system of difference equations 

hm{a--n+m ^ O-n+m) = 0, (7) 

m 

^ ^ [^m(^n— m ^— n+m) ^" Qmis^n+m ^n — m)] — 
m 

which can be solved with an ansatz 

a„ = a+A", 6„ = 6+A", (8) 
a_„ = a_A", 6_„=6_A", 

for n > 0. Quite surprisingly, the resulting homogeneous 
system 



/ h{X) -h{\~^) 
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has rank 2 for any value of the spectral parameter \^ where 
h{X) and ga{^) '■= .9(A) — .^(A"^) are the polynomials 



m-|- 



MA)= 5(A) = 5] 5mA" 



m— — m_ 



m— — m_ 



Hence, there are two linearly independent solutions of (|^) 
(up to an arbitrary common prefactor), namely 



a+(A) = MA-^), b+iX) ^ g{X-'), 
a_(A) = MA), 5_(A)=.g(A), 



(10) 



and 



a+(A) =a_(A) =0, 6+(A) = fe_(A) = 1. (11) 

(i) Let us first consider the case where the solutions 
a±(A),6±(A) are given by (|lO|). The global uniform so- 
lution (for all n G Z) is given by a linear combination of 
TV := m+ + m_ + 1 solutions (p^) 



an = ^ Cma+(A,„)A^, 6„ = ^ c,„6+(A„i)A"j, (12) 



m— 1 

for n > 0, and 

N 



m— 1 



On = ^ Cma_ (Am)A„", ^ C„i6_ (Am)A, 



— n 



m— 1 



(13) 



for ri < 0. iV— tuple of spectral parameters X — 
(Ai . . . Ajv) is an arbitrary subset of a complex unit disk 
Am I < 1 (in order to ensure convergence of T) while the 
coefficients Cm are determined by gluing the solutions 
( |T2| ) and ( p^ ) on m+ + m_ = — 1 sites around n = 0, 
giving a homogeneous system of A^ — 1 linear equations 
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m— 1 



(A;;, - A7„")c„ = 0, n=l...N, 



(14) 



with a general (polynomial) solution 



c™(A) = (-i)™A^r' n(i-^^-^^-) n(^^-^^o (15) 

i<fe j<k 

Logarithmic transfer matrix T{X) is a holomorphic func- 
tion in A, and the coefficients of its Taylor expansion 
around A = also commute with H. After some simple 
series manipulations we easily find an infinite sequence 
of independent integrals of motion, namely the conserved 
charges Qk, k > 0, [H, Qk] = 0, (note that Qo — 2H) 

Qk = ^ [hm{Uk+m + 

k-\-ra ) + 9m{Vk+m k-\-rn )\ 

m— — m — 

(16) 

(ii) In the other case, the solutions a±(A),6±(A) are 
given by (11) and aready solve (Q) globally (so A^ = 1). 
The logarithmic transfer matrix is now rather trivial, 
"^(A) = Sm=i(^ + ^-m)A'", giving the conserved cur- 
rents Ck, fc > 0, [H,Ck] = 0, 

Cfc = Vfe+i + V-k-1 = ■Z'[i(3fc)2] - Z[2(3k)i]. (17) 

Co is the particle current of the associated spinless 
fermion model (via Wigner- Jordan transformation), Ci 
is the energy current, etc. 

It can be easily verified directly that [T(i) (A), T(i) (/!)] = 
0, [7'(i)(A),T(,,)(/i)] EE 0, and [T(,,) (A), T(ij) (/x)] = 
0. Hence all the conservation laws are in involution 
[Qk,Qi] = [Qk,Ci] = [Ck,Ci] = 0. For example, for 
the Ising model in a transversal magnetic field, H = 
JUi + hUo, one recovers well known conservation laws 
Qk = J[Uk+i + Ui-k) + h{Uk + U-k) and Ck which 
required more involved methods in the literature [Q . 
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Conservation laws in non- autonomous case, kicked-Ising 
model: Next we study more general possibly non- 
autonomous quantum spin chains which are propagated 
by members of a unitary Lie group generated by DLA 
6 which in general cannot be written in terms of some 
Hamiltonian H, as exp(— iiJ). For simplicity, we consider 
periodically kicked systems which correspond to time- 
dependent Hamiltonian 

Hit) ^Ho + 6pit)Hi (18) 

where 6p{t) is a periodic delta function of period 1, and 
Hq, Hi G & are the generators — the kinetic energy and 
the potential, respectively. Using the adjoint represen- 
tation of DLA, the (linear) Heisenberg map U^'^ of an 
observable A G & for one time step is factorized as 

A{t + 1) = W^^Ait) = U^'^U^'^Ait) (19) 

where U^'^A — exp(i a.dHp)A — exp{iHp)Aexp{—iHp), is 
the propagation by the kinetic energy and the potential, 
for p — 0, 1, respectively. Transfer matrix is now sought 
by the invariance condition 

U^'^TiX) =TiX) (20) 

in the form (^). The method of solution is analogous 
to (0-^5|) whereas the difference equations for a„ , 6„ are 
now obtained by means of adjoint representation of prop- 
agators which can be derived explicitly by means of eqs. 

and series expansion of exponential function; say if 
generated by Um 

e:X.p{ia adUra)Un = C^Un + S^C/2m-„ + Cs(y„-jn-Vrn-n), 

e:x.p{ia adUm)Vn = c^V„ + s^VL„ - cs([/„+„ -[/,„_„), 

where c = cos(2a), s — sin(2a). 

Here the general procedure cannot be written as ex- 
plicitly as in the autonomous case, so we work out in 
detail an example of kicked Ising (KI) model where the 
kinetic generator is the usual 1-dim. Ising Hamiltonian, 
Hq = JUi = J2j '^'^j'^j+i' ^^"^ kick potential is the 
transversal magnetic field. Hi — hUo = J^j hcr^j ■ Condi- 
tion ( |20| ) results in the system of second order difference 
equations for a„, a_„, 6„, 6_„ which is solved thru the 
ansatz (||) giving the solution (again for any |A| < 1) 

a+(A) = sjCh + cjShX^^, b+{X) = sjSh{X - A"^)/4, 
a^{X) = sjCh + cjShX, 6_(A) = -6+(A) (21) 

where sj = sin(2 J), cj = cos(2J),s/i — sin(2ft,),c/i = 
cos(2/i), and the trivial solution (pi]). In order to ob- 
tain the global solution (a„,5„) we again glue together 
a linear combination of partial solutions for positive and 
negative n at the two points (since the system is of sec- 
ond order), say at n = 0, 1. We obtain the system ( p^ ) 



for the three coeficients c,„, = 3, depending on a triple 
of spectral parameters A = (Ai, A2, A3) with the solution 
dlSl). Collecting the terms with different powers of Am 
in the power series expansion of the logarithmic transfer 
matrix T(A) we obtain two infinite sets of conservation 
laws, namely the charges Qk, k >0, If^Qk — Qk, 

Qk = cjSh{Uk+i + U-k+i) + sjChiUk + U-k) 

- sjSh{Vk+i + V-k+i - Vk-i - V-k-i), (22) 

and currents (0), Ck, k>0, V^Ck = Ck- The conser- 
vation laws of KI model (|2^) are identical to invariants 
( p6| ) of an autonomous IQM system with the Hamilto- 
nian i/Ki = (9o/2 = cjShUi + sjChUQ-sjSh{Vi-V-i)/2. 
Note, however, that the full dynamics arc not identical, 
exp(iadi7Ki) ^V"^. 

Structure of DLA and diagonalization of conserved 
charges: Let us now analyze the structure of DLA & more 
carefully. The current invariants C'k are rather trivial; 
they span a maximal ideal 3 of DLA 6, [6, 3] =0. The 
derived (semi-simple) DLA S' = [6,e] = (adS)°°6 = 
6/3 is spanned by Um ± C^-m and Vm — V-m, for to > 0, 
or in terms of real non-local Fourier transformed basis 
00 

j'i^) = T- E e"^"(C/™ - [/_,„), 

m— — 00 

00 

= ^ E e'™(Kn-V^-™), (23) 

m— — 00 
^ 00 

J'(«;) = -— V e^'""(C/™ + C/-™), 
Stt -'^ 

m—~oo 

for < K < TT, where the commutation relations read 

[JP{k), r{K')] = i5{k - k') J2 eprsJ'in). (24) 

s 

Therefore, derived DLA is isomorphic to an infinite di- 
rect sum & ® 51^2 ■ It has an infinite rank, 
Cartan subalgebra is spanned by a continuous root 
basis J^{k) and Chevalley generators are J^{k) = 
J^(k) ± ij'^^n). Now we construct the vacuum state 
|0) by the condition J~(k)|0) = 0, which is is equiv- 
alent to {Um — U-m — iVm + iV-m)\^) = 0, and also 
crj\%) = 0. Hence the vacuum is the state with all 
spins down. Let us write the Fourier transform of the 
currents and charges as, C(k) = ^„exp(zKTO)Cm and 
Q{i^) = X]jn ^^P(*'*™)Qni ~ '3(~'*); respectively. Note 
that Q{k) — Q{~k) since Qm = Q-m- Using explicit 
form ([l^) we compute 

Q(k) = -87rg(K) • J(k) + 5^(k)C(k), (25) 

q{K) = {hi{K),gi{K),hr(K)), 

where hrin) = Re ft.(exp(iK)), hi{K) = Im ft,(exp(iK)), 
grin) = Re g(exp(iK)), giin) = Im g{exp{iK)). The 
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structure (|2j) is invariant with respect to arbitrary local 
(K-dependent) rotation (non-abelian gauge transforma- 
tion) of the vector field J{k). Particularly interesting is 
the rotation R(k) around axis (?(«;), 



{gi{K), -hi{K),0) 



for an angle <p(k), 



a{K) ■ q{K) = 0, 



(^(k) = arctan 
namely, Rr = (a • f)d + 



where k 



\q\ |ax<j| |axfc| 

(0, 0, 1), which has the property K{K)q{K) — \q{K)\k. The 
unitary transformation of the vector field 

W{k) = R(k)J(k) = ex.p{iJ^dKipa- J)J exp{—iJ^dKipa- J), 

makes the the conserved charges Q{k) proportional to 
the new root basis W3{k) — q{n) ■ J(K)/|g(K)|, namely 



(26) 



Using the same rotation we construct a new vacuum state 
\%)w relative to the field W{k), W-[k)\%)w = 0, namely, 

\^)w = exp J dKip(K)d{K) ■ J(k)^ |0). 

Let us now discretizc the momentum k to L bins what 
corresponds to (but is not identical to) a finite chain of 
L spins, and define 



TTfc/L 



7r(fc-l)/L 



dKW^in), l<k<L. 



Then we have [Wj^,W[] = i5kiT,,eprsW^. The eigen- 
states of conserved charges Q{k) (and of root basis 
since [Q{k),W^{k')] = 0) can be labeled by L binary 
quantum numbers G {0, 1}, k — 1 . . . L, and are con- 
structed by means of creation operators 



Ck = l 



ic,)= n ^^m 



(27) 



1<A;<L 



with Wj^\ch) — [ci — l/2)|cfc). Hence, for smooth \q{n)\ 
and large L, all the charges ( p6|) are diagonal in the 
eigenbasis (p7| ) . Of course, the eigenvalues are finite (for 
the infinite system L = oo) only for the charge densi- 
ties Q'j^ = limi^oo(l/-^)'9m|L and not for the extensive 
charges Qm- The eigenvalues of invariant densities are 
computed by taking the limit L ^ oo and the inverse 
cosine transform 
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Q'„^\c{k)) — / dK.'cos{n' ■m)\q(n' 

71" Jo 



|c(^)). 



Label c{k) is an arbitrary (measurable) index function 
c : [0,7r) ^ {0,1}, (for finite L, c{k) = Cfc,7r(fc - 1)/L < 
K < nk/L) and \c{k)) is the corresponding eigenstate 
which should be properly defined by some limiting pro- 
cedure L — > oo of (|2^) . It seems that the spectrum of Q'^ 
is purely continuous. On the other hand, for the currents 
we have Cm\c{K)) = 0, since [C{k),Wi] = 0. Having such 
a transparent structure ( p4|^7| ) it should be an easy task 
to compute physically interesting correlation functions. 

In this Letter we have introduced an infinitely dimen- 
sional space of completely integrable quantum many- 
body systems (spin-chains or chains of spinless fermions), 
the so called Dynamical Lie Algebra, as opposed to 
few-parameter families of completely integrable quan- 
tum many body systems known so far in the literature. 
The model is an infinitely dimensional extension of Ising 
model in transversal field (equivalent to XY-modcl Q| 
and to 1-dimensional free fermion theory). For every el- 
ement of the algebra being interpreted as a Hamiltonian, 
or any propagator from the associated unitary Lie group 
being generated by a finite number of elements of the al- 
gebra (e.g. Ising model periodically kicked by transversal 
magnetic field), we construct two infinite sets of quan- 
tum invariants of motion, the conserved charges and the 
conserved currents. Is is shown heuristically how to di- 
agonalize these conservation laws. Explicit expressions 
of the conserved charges are quite simple (much simpler 
than in general Heisenberg (XYZ) or Hubbard model 
for example) though nontrivial. 
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